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We solve the Firing Squad Synchronization Problem (FSSP), for P systems based on digraphs with 
simplex channels, where communication is restricted by the direction of structural arcs. Previous 
work on FSSP for P systems focused exclusively on P systems with duplex channels, where commu- 
nication between parents and children is bidirectional. Our P solution, the first for simplex channels, 
requires cell IDs, strongly connected digraphs and some awareness of the local topology (such as 
each cell's outdegree) — we argue that these requirements are necessary. Compared to the known 
solutions for cellular automata, our solution is substantially simpler and faster. 

Keywords: P systems, digraphs, strongly connected, simplex channels, firing squad synchronization, 
cellular automata. 

1 Introduction 

The Firing Squad Synchronization Problem (FSSP), originally proposed by Myhill in 1957 lfl4l . is one 
of the best studied problems for cellular automata. Essentially, the problem involves programming a 
network of cellular automata, so that, after the firing order is given by the general, after some finite time, 
all the cells in the system enter a designated firing state, simultaneously and for the first time. 

Versions of FSSP have been proposed and studied for variety of network structures, from simple lin- 
ear graphs to rings, trees, or general connected graphs; see, for example, BI^ [T0l[TT1[T^[T31[T7ll23ll25l 
[26). However, most of these versions, require bidirectional communication (i.e. duplex channels): only 
a few notable exceptions consider the more restricted unidirectional communication (i.e. simplex chan- 
nels), starting with Kobayashi lfl2l . Later, Even, Litman and Winkler [9] proposed improved solutions 
for arbitrary undirectional strongly-connected digraph, working in 0(N 2 ) steps, where N is the digraph 
size (number of cells). Ostrovsky and Wilkerson [18] improved this further, to a solution which runs in 
0(ND) steps, where D is the digraph diameter (typically smaller than AO — this still seems to be the best 
solution available. 

Several FSSP solutions have recently been studied in the framework of P systems, although with 
somehow different formulations, stemming from their different computing capabilities. P solutions were 
proposed: for trees, by Bernardini et al. [3] and Alhazov et al. [1]; and for arbitrary connected graphs, 
by Dinneen at al. Ed [8). All these P solutions require duplex channels and follow the typical pattern 
of a wave algorithm E4l . using three phases: 
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1. a first broadcast — which follows all shortest paths from the general and builds a virtual BFS tree 
(or dag); 

2. a convergecast — which helps determine the general's eccentricity; 

3. a second broadcast — which carries the actual firing command (with a countdown counter). 

The best P solutions need e g +k steps for each of the three phases, for a total of 3e g + k steps, where e g 
is the general's eccentricity (height for trees, if the general is at the root). 

Obviously, while the two broadcasts, of phases (1) and (3), would also work with simplex channels, 
duplex channels are essential for the convergecast of phase (2), where children need to talk back to their 
parents. At first sight, the convergecast seems impossible for simplex channels. However, children can 
still talk back to their parents, if the digraph is strongly-connected, albeit on a typically longer path. 
Moreover, if messages cannot be confused, all children can send messages to their parents, in parallel 
(overlapping in time without problems), achieving this way a virtual convergecast. 

Based on ideas from Ostrovsky and Wilkerson lfl8l . we propose & first FSSP solution for P systems 
with simplex channels, based on arbitrary strongly -connected digraphs. Our solution runs in 0(e g D) 
steps, specifically: (1) e g steps for the first broadcast; (2) 0{e g D) steps for the virtual convergecast 
(maximally parallelized); (3) e g steps for the second broadcast. Thus, in terms of execution time, our 
P solution compares favourably with the best known solution for cellular automata. Taking into account 
the different problem constraints and different computing capabilities of P systems vs. cellular automata, 
we were expecting a simpler and faster solution, but not necessarily such a substantial speed improve- 
ment. However, for the reasons mentioned, any performance comparison must be viewed with a grain of 
salt. 

The actual design seems challenging and requires careful selection of the most adequate ingredients, 
some of which are available in cellular automata, but not typically available in P systems. Specifically, 
we argue that the P solution requires: (1) reified cell IDs; (2) reified local network information, such 
as the number of outgoing arcs (this information is available in cellular automata); and (3) high-level 
generic rules, if we want to have a fixed rule set, independent of the actual network size. 

2 Preliminaries 

We assume that the reader is familiar with the basic terminology and notations, such as relations, graphs, 
nodes (vertices), edges, directed graphs (digraphs), directed acyclic graphs (dags), arcs, alphabets, strings 
and multisets. 

A P system is a parallel and distributed computational model, inspired by the structure and interac- 
tions of cell membranes. This model was introduced by Paun in 1998-2000 [19]. An in-depth overview 
of this model can be found in Paun et al. |j22j. 

In this paper, we consider an ad-hoc definition of P systems, based on our definition of simple P mod- 
ule 0, which extends earlier versions of tissue and neural P systems iPTBIEOll . However, here we inten- 
tionally restrict rule transfer mode to broadcast to all children, |v 

Definition 1 A P system of order n with simplex channels is a system IT = (0,K,8), where: 

1. O is a finite alphabet of elementary symbols; strings over O are interpreted as multisets; 

2. K = {(7i,(72,...,(7 n } is a finite set of cells; where each cell is a system Oi = (Qi,Ri), with Qi a 
finite set of states and Ri a finite set of rewriting rules over O, further detailed below. 
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3. 5 is an irreflexive binary relation on K, which represents a set of structural arcs between cells, 
with unidirectional communication capabilities, strictly from parents to children. 

4. Each Ri is a finite linearly ordered set of multiset rewriting rules with promoters, of the form: 
S x — > a S' x 1 (y)i v ■■■ \z> where S,S' £ Qi, x,x' ,y £ O*, z£ O* is the promoter, a G {min,max} is a 
rewriting operator and |v is a transfer operator, here restricted to sendy messages from a parent 
to all its children. 

As usually, each cell, a, G K, starts from an initial configuration (5,o,w ; -o), where £ Qi is its initial 
state and w,o 6 O* is its initial content. A cell evolves by applying one or more rules, which can change 
its current configuration, i.e. its current state and current content, and send symbols to its children. 

The application of a rule transforms the current state S to the target state S', rewrites multiset x as x 1 
and sends multiset y by replication to all its children. Note that, multisets x' and y will not be visible to 
further rules in this same step, but they will become visible after no more rules are applicable, i.e. they 
will be available since next step only. Promoters are symbols which enable rules, but are not consumed 
by the rules' application. 

When an applicable rule is applied, its rewriting operator a indicates how many times it is actually 
applied: once, if (X = min; or as many times as possible, if a = max. 

As used here, rules have priorities and are applied in weak priority order ll2TTl . with special attention 
to target state compatibility: (1) higher priority applicable rules are applied before lower priority appli- 
cable rules, and (2) a lower priority applicable rule is applied only if it indicates the same target state as 
the previously assigned rules (if any). 

All cells evolve synchronously in one global step. An evolution of a P system is a sequence of steps, 
where each cell starts from its initial configuration. An execution halts if no cell can evolve. 

2.1 Further P systems extensions 

We let each cell, a,, start with its own unique cell ID symbol, l;. We thus reify the conceptual cell index, i, 
into an internal symbol, which is accessible to the rules, exclusively as an immutable promoter iPToll . 

We enhance our vocabulary by recursive composition of elementary symbols from O into a simple 
form of complex symbols Ifl6l . Such complex symbols can be viewed as complex molecules, consisting 
of elementary atoms or other molecules. 

Further, complex symbols let us process our multisets with high-level generic rules, using free vari- 
able matching. To explain these additional ingredients, consider this hypothetical rule (which uses an 
additional transfer mode, targeted to a specific child, not considered in the definition used here): 

S a nj hnin.min $ b (Cj)^j \lj. 

This is a generic rule, which uses an extended rewriting mode, with complex symbols, c, and nj, 
where i and j are free variables. In fact, c, and nj are just shorthands for tuples (c,i) and (n,j), or, 
equivalently, for compound terms c[i) and n(j). If needed, we can build more complex symbols by 
recursive composition; e.g., we could have complex symbols such as d(e,i,f(j)). Generally, a free 
variable could match anything, including another complex symbol. However, in this rule, i and j are 
constrained to match cell ID indices only: 

1. i — because it also appears as the cell ID of the current cell, if, 

2. j — because it also indicates the target of the transfer mode, 
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A generic rule is identified by using an extended version of the "classical" rewriting mode, in fact, by 
a combined instantiation and rewriting mode. Our sample rule uses the extended mode min.min, where 
the two min operators have distinct semantics: the first min operator is new and describes the generic 
instantiation; the second min is the classical operator, which describes the rule application. Briefly: 

1. according to the first min, this rule is instantiated once, for one of the existing nj symbols (if any), 
while promoter, l,-, constrains i to the cell ID index of the current cell, Of, 

2. according to the second min, the instantiated rule is applicable once, i.e. if applied, it consumes 
one a and one tij, produces one b and sends one c; to child Oj (if this exists). 

As a further example, consider the scenario in which the current cell, 0\ , contains the multiset n^n^n^. 
Here, our sample generic rule instantiates (randomly) one of the following two lower-level rules, which 
is then applied in the classical way (in the min rewriting mode): 

San 2 -WS' b (ci); 2 . 

S'b( Cl ) h . 

We consider four basic combinations of the instantiation and rewriting modes, min.min, min.max, 
max. min, max. max; their semantics is: 

• min.min indicates that the generic rule is (randomly) instantiated once, if possible, and the instan- 
tiated rule is applied once, if possible. 

• min.max indicates that the generic rule is (randomly) instantiated once, if possible, and the instan- 
tiated rule is applied as many times as possible. 

• max.min indicates that the generic rule is instantiated as many times as possible, without superflu- 
ous instances (i.e. without duplicates or instances which are not applicable), and each one of the 
instantiated rules is applied once, if possible. 

• max.max indicates that the generic rule is instantiated as many times as possible, without superflu- 
ous instances (i.e. without duplicates or instances which are not applicable), and each one of the 
instantiated rules is applied as many times as possible. 

All instantiations are ephemeral, created when rules are tested for applicability and disappearing at 
the end of the step. 

3 FSSP problem for P systems with simplex channels 

We are required to find: 

1. an alphabet O; 

2. a cell prototype a = (Q,R), where 

(a) R is a set of rules over O; 

(b) Q contains two distinguished states: 

• So: a quiescent state, i.e. if a is in state So and empty, then there are no applicable rules; 

• Sf: afinal state, i.e. if a is in state Sf, then there are no applicable rules. 

such that, given: 

1. any finite set of a copies, K = {o\ , 02, ■ ■ . , o n }, a, = a; 



60 



2. connected via any strongly-connected digraph <5; 
the P system IT = (0,K, S) with simplex channels will evolve according to the following specification: 

1. all cells start from quiescent state So: S;0 = So; 

2. except a distinguished cell a g , called the general, all cells start with a restricted initial content, con- 
taining, the reified cell ID and a reified count of the cell's outdegree: w i0 C {ij,c outdegree ^''}, V/ / 
g\ 

3. the evolution terminates and, during its last step: all cells enter state Sf simultaneously and for the 
first time. 

Remark 1 Our formulation has different constraints than the original problem for cellular automata. 
In the cellular automata formulation, there is a given fixed bound on the number of input and output 
connections of each cell (bounded indegree and outdegree). In our formulation, there is no bound on the 
number of input and output connections that a cell may have. However, this is compensated by the fact 
that in our formulation there are no size bounds on messages or cells ' internal memory. These trade-offs, 
as well as different computing capabilities, suggests that performance comparisons must be viewed with 
a grain of salt. 

Remark 2 We argue that both the reified cell ID and the reified children count, or equivalent informa- 
tion, are necessary, definitely for our approach, and, likely, for any other approach. Note that children 
counts are implicit in the cellular automata version, where unconnected channels (out of the fixed sized 
pool) can be detected. 

Remark 3 There is no constraint on the cells' final contents. However, if needed, any left-over garbage 
could be collected in one extra step. 

Remark 4 Practically, we are only required to design the rule set, R, because this implies the alpha- 
bet, O, and the state set, Q. 

Remark 5 Note the rule set, R, must be fixed and applicable to any structural digraph. This is a strong 
requirement: we require a rule set which is independent of the size and structure of the actual system. 

4 FSSP solution for P systems with simplex channels 

Our solution runs in three phases (conceptually similar to the duplex case): 

1. First phase: & first broadcast from the general. This phase builds the virtual-dag (the virtual BFS 
dag) and, for each cell: (a) records its virtual-dag-parent(s); and (b) successively computes its 
depth attribute, which represents this cell's depth level in the virtual-dag (the same as this cell's 
digraph distance from the general). A first phase broadcast message is a complex symbol, x^j, 
where i is the sender's ID and k is the next depth level (a,'s own depth plus one). 

2. Second phase: a virtual convergecast from the virtual-dag leaves. For each cell, this phase succes- 
sively computes the max-depth attribute, which represents the maximum depth over all descendant 
cells in the virtual-dag. In the end, the general's max-depth is its eccentricity. 

This virtual convergecast simulates impossible direct virtual-dag-child to virtual-dag-parent mes- 
sages, by broadcasting them over the digraph (using ad-hoc BFS dags). A convergecast message 
is a complex symbol, a/,;,*, where i is the sender's ID, j is its virtual-dag-parent ID and k is a,'s 
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own max-depth. In addition to virtual-dag leaves, a pseudo-convergecast message is also sent by 
a digraph cell to its digraph parents with larger depth attributes (if any). Because each message is 
uniquely identified by both it sender and its destination, any number of such convergecasts can run 
in parallel, without creating confusions. 

Note that a cell needs to know its digraph outdegree — to detect when it receives its last outstand- 
ing convergecast message and to start its own convergecast. However, a cell does not know the 
identities of its digraph children, or how long a message from any one of them will take to reach 
it. 



3. Third and last phase: a second broadcast from the general, with a countdown to firing. A last 
phase broadcast message is a complex symbol, fa, where k is the next countdown counter (a,'s 
own countdown minus one). 



As possible extensions, not discussed here, we can extract from our solution a more general subpro- 
gram, to send any message from any cell to any other cell, which can run in parallel, without creating 
confusions. Also, we can consolidate the routing information, to speedup future messages with the same 
destination; however, this feature is not required here (each convergecast is performed exactly once). 

Figures [T]-[8] show bird's eye views of the evolution of II: a sample P system, with simplex channels, 
based on a strongly connected digraph. 








*8l- 1- 




*9l- I" 




o 







CT 10 


-1- 



Figure 1: IT: Initial configuration of a sample P system with simplex channels, based on a strongly 
connected digraph. Normal arrows are arcs in the virtual dag, created by a BFS broadcast, started from the 
general, <3\. The remaining digraph arcs are dotted arrows. Each cell knows its digraph outdegree, indicated by 
small circles at outgoing arrows' tails. Each cell blob shows three attributes, in order: (1) its cell ID; (2) its depth 
attribute (computed by the first broadcast); and (3) its max-depth attribute (computed by the virtual convergecast). 
At this stage, the depth and max-depth attributes are still indeterminate and all small circles are white, indicating 
still outstanding convergast messages. 
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4.1 Rule set R 



0. Rules in state So- 

1. Sq Cl hnin.min 

Si gn m (xu); v \U 

2. Sq X k j ^max.min 

Si h Pj (^+i,i); v \ l i 

3. Sq X k J ^max.max Si Pj \li 

1. Rules in state Si: 

1. S() h — ^max.min Si fl k m k 

2. S() ^max.max Si 

3. Si XfcJ ^max.min Si y j.Q 

4. Si Xfcj ^max.max Si 
J. Si S k S'max.min 

Max w m k 

6- Si yjfi ^max.min 

51 Vj.i (a ; -,a); v \M 

7. Si Clj,k,l ~ ^max.max Si \Vj£ 

8. Si Oijjk ~ ^max.min S] V,j 5/; |lj 
9- Si O/^/ — )- max .mm 

51 V;,* (a^,;); v 

2. Rules in state S2: 

1. S2 ? — ^min S5 (7); v 

2. S 2 *jt -^max.min MdX W m k 

3. S 2 yj,k — ^max.min 

52 v,y {ciju)u \h 

4. S2 a,- 1,1 r max.max 

5. S2 a; it nnax.min 

6. S 2 a j,k,l ^max.min 

5 2 v,-,* (aj^i)i v 

3. Rules in state S3: 

1. S3 cc 

'min 

Si ce 

2. S 3 c -> min S 4 eft 



3. Rules in state S4: 

1. S4 ^max.min 

S2 1 f k (0;v Is 6 m <t 

2. S 4 ^max.min 

4. Rules in state Max (needs re- 
finement): 

1. Maxm k -> max S 3 \m k+ i 



5. Rules in state S5 : 


1. 


S5 v u 


^max.max Sf \/q 


2. 


S5 w 


^max.max Sf \fo 


3. 


S5 


~~ ^ max. max Sf \fo 


4. 


S 5 


^max.max Sf \/q 


5. 


S 5 e - 


^min.max Sf \fo 


6. 


s 5 6- 


"^min.max Sf \fo 


7. 


s 5 -> 


min.min Sf \fo 


8. 


S5 Vfc,/ 


^max.max S5 \f k 


9. 


S5 At 


^max.max S5 \f k 


10. 


S5 


^ max. max S5 \f k 


11. 


S 5 


^max.max S5 \f k 


12. 


S 5 e - 


^min.max S5 \ f k 


13. 


s 5 6- 


"^min.max S5 \f k 


14. 


$>/*- 


'max. mill 






S5 fk-l (fk-l)i v 


15. 


s 5 A- 


~ ^max.max S5 


16. 


s 5 ?- 


^max S5 



4.2 Alphabet O, elementary and complex symbols 

• if. reified cell ID; 

• a: starts the process from the cell which next assume the general role; 

• g: marks the general; 

• xj£. complex symbol broadcasted in the first phase; 
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• n^. indicates the depth attribute, k; 

• mf. indicates the max-depth attribute, k; 

• Ik', auxiliary symbol used to compute depth attribute, 

• pk'- pointer to a parent cell; 

• yi£. initiates the sending of message k to Of, 

• ciijy. complex symbol broadcasted in the convergecast phase, sent from Oj to a, and carrying 
payload k, representing a/s max-depth (if known, otherwise it is a pseudo-convergcast sent by a 
cell with lower depth); 

• Vjj: records the passage of a message from O) to Oj; 

• auxiliary symbol used to compute max-depth attribute, m^, via maximum; 

• c: used to count the children which have not yet sent their convergecast messages (initially cell's 
outdegree); 

• e: used to count the children which have already sent their convergecast messages (initially zero); 

• b : marks a cell which has performed its convergecast; 

• t: auxiliary symbol used in the countdown to firing; 

• fu'- complex symbol broadcasted in the last phase, carrying the countdown to firing; 

4.3 Brief description 

1. State So: initiates the first broadcast and computes depth attributes. 

2. State S\ : continues the first broadcast and builds the virtual-dag. 

3. State S2- performs the actual convergecast. 

4. State S3: decides if a non-general cell is ready to start its convergecast, i.e. if it has received its all 
outstanding convergecast messages, from all its digraph children. 

5. State S4: decides if the general is ready to start the second broadcast (countdown to firing), i.e. if 
it has received convergecast messages from all its children. 

6. State Max: prepares the convergecast, by determining the max-depth attribute. 

7. State S5: last broadcast, counts down to firing and erases unnecessary symbols. 

5 Assessment 

Theorem 1 The synchronization time of the FSSP solution for digraph-based P systems with simplex 
channels is e g + e g D + e g , i.e. bounded by 0(e g D). 

Theorem 2 The digraph must be strongly connected, otherwise, there is no solution. 
Theorem 3 Cells must know the number of their outgoing arcs, otherwise, there is no solution. 
Theorem 4 Cell IDs must be reified, otherwise, there is no solution. 
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6 Experimental results 



Besides our earlier example, we have empirically validated our solution in several test scenarios, with 
digraphs of different shapes and sizes, for example: 

• A simple ring (Figure [9]). 

• A main ring linking a series of smaller rings of size two (Figure 10). 



• A main ring linking a series of smaller rings of size three (Figure 11). 

• A main ring linking a series of smaller rings of increasing size (Figure [T2|). 

• A set of 1 1 random directed graphs, with up to 70 nodes each, generated by the standard networkx 
package. 

The results support our claim for correctness and performance. 

7 Conclusions 

In this paper, we explicitly presented a first solution to the FSSP for synchronous digraph-based P sys- 
tems with simplex channels. Our design suggests, but does not need, ways to consolidate routing infor- 
mation in such systems — this can be a topic for further study. 

Our solution runs in O(eD) steps and compares favourably, i.e. it is simpler and faster, than the 
best known solution for cellular automata lTT8ll . which runs in 0(ND) steps. Taking into account the 
different problem constraints and different computing capabilities of P systems vs cellular automata, we 
were expecting a simpler and a faster solution, but not necessarily such a substantial speed improvement. 
However, as noted before, any performance comparison must be viewed with a grain of salt. 

Our solution used a fixed size high-level rule set, independent of the number of cells in the actual 
system and of its structure. This supports the case for reified cell IDs, complex symbols and generic 
rules and suggests that such ingredients could be useful or even essential in any distributed or just large 
system. 
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Figure 2: IT: route from Oj to 02, indicated via thick arrows. 

joi- I 




°"2|-|- CT 3|-|- g 4l - I 



Oo-gl-l- <T 6 |-|- CT7I-I- 








=■"10 


-1- 



Figure 3: IT, step 1: The general, G\, starts the first phase, by broadcasting the complex symbol ai t 
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Figure 4: IT, step 9: Cell <7io has just learned that it is a virtual-dag leaf, by receiving cell's 04's pseudo- 
convergecast message, the complex symbol aio,4,-. Cell O\o starts now its virtual convergecast, by 
broadcasting the complex symbol 09,10,4, towards its virtual-dag-parent, G9. At this stage, each cell knows 
its virtual-dag-parent(s) and its own depth attribute. All other cells have already initiated their virtual convergecasts. 
All leaves, including O10, and some other cells already know their max-depth attribute: exactly, if they have 
received all their outstanding convergecast messages, or a lower bound, otherwise. Message 09.10,4 will take three 
more steps: via paths (7] 0.04.06. 09 and path 010-04. 07.O9. Three other, earlier stared, virtual convergecasts run in 
parallel: 02 to 01, 03 to Oi, 0% to a^. Smaller arrows near structural arcs indicate virtual convergecasts. 
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Figure 5: IT, step 12: After receiving its single expected convergecast message, cell 09 starts its own 
virtual convergecasts towards its parents, at and 07, by broadcasting complex symbols, 06.9.4 and 07,9,4, 
respectively. Each of these messages will take three more steps, via path O9.O10.O4.O6 and via path 09. 010-04. 07, 
respectively. These two convergecasts run in parallel with another virtual convergecast: Os to 06- 
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Figure 6: IT, step 18: The convergecast phase is almost completed. Cell 04 starts its convergecast 
towards its parent, Oj, by broadcasting the complex symbol a 1.4,4. This convergecast will take four more 
steps, via paths a4.GQ.a3.a5.a1 and path c74.C76-C78.C75.C7i. This is cell C7i's last oustanding convergecast. 
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Figure 7: IT, step 22: General Oi starts the last phase, by broadcasting the complex symbol fa, carrying 
the countdown to firing. Small dots above a cell indicate the cell's countdown counter, which is also broadcasted 
to all its digraph children. 
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Figure 8: IT, step 25: The last phase is almost complete. Cell 09 forwards the complex symbol fo to 
ClO. This the last step before firing (not illustrated here). 




(a) Sample network, N = 6. 
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(b) Results. 
Figure 9: Ring networks. 
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(a) Sample network, N = 10. 
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(b) Results. 

Figure 10: Ring networks of size 2 rings. 




(a) Sample network, Af = 18. 
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Figure 11: Ring networks of size 3 rings. 




(a) Sample network, Af = 14. 



Figure 12: Ring networks of increasing size rings. 
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